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Abstract 

We study all tree-level neutralino annihilation processes concerning energy-dependence at the 
mSUGRA benchmark point SPSla using tree-level calculators MadGraph and CompHEP as well 
as FormCalc. We further investigate the total tree-level neutralino annihilation cross section by 
scanning over the mSUGRA parameters and identifying the dominant channels in every region. 
The results for the cross sections are presented in the m 1 / 2 - , ma planes for different values of tan/3 
and A . We identify as well the resonances due to the exchange of the heavy Higgs bosons A 
and H in the s-channel of many neutralino annihilation processes, where the total cross section 
increases significantly if these Higgs bosons can be produced close to on-shcll. Finally, we provide 
a few comments on the neutralino annihilation processes at loop-level. 
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1 Introduction 



In recent years great efforts have been performed in the development of theories beyond the standard 
model of particle physics. One of the most promising theories is supersymmetry (SUSY). In super- 
symmetric models many new particles are introduced and some of them are expected to be found at 
the upcoming collider LHC at CERN. 

Supersymmetry is able to solve the most severe problems of the standard model in a quite natural 
way. In addition it may supply a promising candidate for dark matter under the assumption of an 
additional discrete symmetry, i.e. R-parity conservation. Under this symmetry SUSY-particles can 
only be created and annihilated in pairs, so that supersymmetric decay chains lead to a stable lightest 
supersymmetric particle (LSP), which is the candidate for dark matter. In most models this LSP is 
the lightest neutralino. Assuming this particle to be the weakly interacting particle the dark matter 
of the universe is made of, the results from measurements of the dark matter relic density from the 
cosmological observations by WMAP pQ can be used to put constraints on supersymmetric models 
[21 [3]. For these purposes it is important to know how the neutralinos may annihilate themselves. 

In this report we study the cross section of neutralino annihilation processes at tree-level in the 
model of minimal supergravity (mSUGRA) , which we review briefly in section [2J After some remarks 
on the neutralino annihilation cross sections and the connection to cosmology in section [3j we focus 
on annihilation processes at tree-level in section Here we discuss first the possible final states and 
Feynman diagrams (|4.ip . Then a brief review on the tree-level calculators CompHEP [4] and Mad- 
Graph [5 J is given (|4.2p . We study the different tree-level neutralino annihilation processes at the 
mSUGRA benchmark point SPSla in section 14.31 concerning energy dependence and channel domi- 
nance. The mSUGRA parameter dependence of the total annihilation cross section at tree-level is 
analyzed in section 14.41 Here we point out the main contributing channels for specific regions in the 
parameter space. In Appendix A we study in detail the contributions and interferences of the different 
diagrams in the process XX ~^ Z°Z°, trying to get a better understanding of the energy dependence 
plots presented in section 14.3.21 We summarize our results in section [5] and give an outlook on future 
investigations including loop-level corrections. 

2 Definition of the supersymmetric model 

This study on neutralino annihilation processes is evaluated and interpreted in the minimal super- 
gravity (mSUGRA) model. This model is defined by only 4 independent parameters plus a sign (in 
addition to the 18 standard model parameters): 

m , m 1/2 , A , tan/?, sign(/i) (1) 

where mo is the soft symmetry breaking scalar mass (universal for all flavors, at the scale Mqut 
of Grand Unification), m\/2 the universal supersymmetry breaking gaugino mass, Aq the universal 
supersymmetry breaking trilinear scalar interactions, tan/3 the ratio of the two Higgs doublet vacuum 
expectation values and sign(^) the sign of the supersymmetric higgsino mass parameter. 

In this model, R-Parity is assumed to be conserved, i.e. there is a discrete symmetry R = (— l) 3 ( B_L )+ 25, 
of the theory, where B and L are baryon and lepton number operators and S the spin. This imposes, 
that R = 1 for standard model particles and R = — 1 for their superpartners, which leads consequently 
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to a stable lightest supersymmetric particle (LSP), since the decay of SUSY-particles to standard model 
particles is forbidden. In mSUGRA, the LSP is most likely the lightest neutralino x?> in the following 
sometimes only denoted by x- 



3 Neutralino annihilation cross sections 

The annihilation of neutralino pairs plays a central role in neutralino cosmology. Its cross section is 
needed for calculations of the cosmological neutralino relic abundance, the flux of energetic neutrinos 
from neutralino annihilation in the Sun and Earth and fluxes of anomalous cosmic rays produced by 
neutralino annihilation in the galactic halo. For these purposes, it is generally sufficient to expand! the 
annihilation cross section a a in the nonrelativistic limit (v — > 0, where v is the neutralino-neutralino 
relative velocity), 

a A v = a + bv 2 + C(v A ), (2) 

where a is the s-wave contribution at zero relative velocity and b contains contributions from both 
the s- and p-waves. The relative velocity of neutralinos in the galactic halo, sun and earth is of order 
0(1CP 3 ), so only the a term in Eq. ([2]) is needed for calculations involving relic neutralinos. When 
neutralino interactions freeze out in the early Universe, their relative velocities are approximately 
v ~ so both the a and b terms are generally needed for relic-abundance calculations. For more 
accurate calculations the term of order 0(v A ) is also included, but generally, one calculates the thermal 
average 

where pi is the three- momentum of the incoming particle i, fi is its equilibrium distribution function 
and Vij is the relative velocity of the interacting particles [7] . 



4 Neutralino annihilation processes at tree-level 
4.1 Neutralino annihilation processes 

There are numerous final states into which the neutralinos can annihilate. The most important of 
these are those which appear at lowest order in pertubation theory, i.e. the final states which occur at 
tree level. Specifically, these are fermion-antifermion pairs // (where / are the standard model leptons 
and quarks), W + W~ , Z°Z°, W + H~, W~ H + , Z°A°, Z°H°, Z°h°, H + H~ and all six combinations 
of ^4°, h and H°. There are various Feynman diagrams contributing to every channel and the cross- 
sections for all these tree- level processes have been calculated and are summarized in [8]. Moreover, 
automatic calculators such as CompHEP [3], MadGraph |5j and FormCalc [9] are able to calculate 
these cross-sections since the MSSM is implemented in these programs. In the following sections we 
will take a closer look at the different final states which are classified as fermionic final states, weak 
gauge-boson final states and final states containing Higgs bosons. 

2 This approximation is valid unless the neutralino annihilation takes place near a pole in the cross section, as shown 
in0. 
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4.1.1 Fermionic final states 



The annihilation of neutralinos to a fermion-antifermion pair reveals several important features. First, 
with the expectation that the neutralino mass m x is of order or greater than 10 GeV, light fermion 
final states will always be an open annihilation channel. For many interesting neutralino masses, other 
channels will be closed or suppressed, so that fermionic final states are often the only open channels. 
Second, there is a suppression of the cross sections to fermionic final states in the limit of zero relative 
velocity due to helicity constraints: since neutralinos are Majorana particles, i.e. they are their own 
anti-particles, their spins have to be opposite in a relative s-wave due to Fermi statistics. Therefore, 
the fermions in the final state must have opposite spins as well which implies that the amplitude for 
this process carries a factor of the fermion mass mj due to the helicity flip. Thus, the cross-section 
is of order mj/m x . So the suppression of light fermionic final states is due purely to the fact, that 
these fermions are relatively light compared to the energy scale 2m x of the process. Of course there 
is no suppression of the top-quark final states if it is open. In fact, in models where the neutralino 
is heavier than the top-quark, this annihilation channel is expected to be the dominant one since the 
cross section depends on the square of the fermion mass. 

* ? \ \ , ^ 

x? / x!^^ x?^ ^/ 



\^ ^/ x!^ ^ 



Figure 1: Diagrams contributing to neutralino annihilation into fermions. 

The diagrams for the neutralino-neutralino annihilation to the // final state are shown in Fig. [TJ The 
cross section contains contributions from t- and u-channel exchange of all sfermions / which can couple 
to / and from s-channel exchange of Z° and all three neutral Higgs bosons. In general, there are six 
sfermions which share the same charge and mix among themselves. All these sfermions can therefore 
couple to / in the t- and u-channels, but in our case it is sufficient only to consider the exchange of 
the superpartner of the right- and left-handed final state fermions. Since there are s-channel diagrams 
exchanging massive particles we expect resonance effects if the energy is sufficient to produce these 
particles on-shell. However, we will approximate all fermions except those of the third generation as 
massless. Since the Higgs particles couple to the mass, s-channel exchange of these will only happen 
for the third generation. 

4.1.2 Weak gauge-boson final states 

Annihilation to weak gauge-bosons will occur if the neutralino mass is high enough for such channels 
to be open (m x > mw) ° r if the colliding energy is high enough. There is no s-wave suppression 
mechanism for these processes and thus they can be very important in the limit v — > if the neutralino 
is heavy enough to produce weak gauge-bosons in the final state. 

The diagrams for the neutralino-neutralino annihilation to weak gauge-boson final states W + W~ and 
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Figure 2: Diagrams of the process XX ~* W + W . 




Figure 3: Diagrams of the process XX ~^ Z°Z°. 

Z°Z° are given in Fig. [2] and Fig. respectively. Annihilation to W + W pairs occurs by s-channel 
exchange of Z°, h° and H° bosons and by t- and u-channel exchange of charginos. Neutralinos 
annihilate to Z°Z° pairs via t- and u-channel exchange of neutralinos and by s-channel exchange of 
H° and h° bosons. In the v — ► limit the cross section is determined solely by the t- and u- channel, 
i.e. the exchange of charginos and neutralinos for and Z°Z° respectively. 

4.1.3 Final states containing Higgs bosons 

The annihilation of neutralino pairs to Higgs bosons may also be important when such channels are 
open. However, there are some cases where the s-wave amplitude vanishes because of the inability of 
the Higgs bosons to produce the appropriate final state quantum number CP = —1. These cases are 
XX ~^ h°h°, XX ~^ H°H°, XX ~^ A°A°, XX ~^ H + H~ and the mixed final state XX ~^ Z°A°. The 
processes which have non-vanishing s-wave contributions are XX ~^ h°A°, XX ~^ H°A°, XX ~^ Z°hP, 
XX -» Z°H° and X X W ± H^. 



\i - 



Xl 



x!^ 



Figure 4: Diagrams for the processes XX ~~ * Z°h° and 



XX 



Z°H°. 



\i - 
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Figure 5: Diagrams for the process XX ~^ Z°A°. 
The diagrams for the processes XX ~> Z°h° and XX ~> Z°H° and the process XX ~> 

Z°A° are shown 

in Fig. 0] and [5j respectively. All these processes receive contributions from t- and u-channel exchange 
of all four neutralinos. In the processes containing a neutral CP-even Higgs boson, h°, H°, in the final 
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state s-channel exchange of Z° and A contribute, while the process XX ~^ Z°A° contains contributions 
from s-channel exchange of h° and H°. 

The diagrams of the process XX ~> W + H~ (and its charge conjugate) are given in Fig. [6l It receives 
contributions from chargino exchange in t- and u-channel and all three neutral Higgs bosons in the 
s-channel. 



-W+ 



:> 



h, H 



-W+ 



> 



Figure 6: Diagrams for the process XX ~* W + H (the charge-conjugate process is similar). 

Annihilation to the two Higgs-boson final states hPhP, H°H°, h°H° and A A occurs via the exchange 
of neutralinos in the t- and u-channel and the exchange of h and H° in the s-channel, as shown in 

Kg. m 



.■■■A 



Figure 7: Diagrams for neutralino-neutralino annihilation to the h°h°, H°H°, h°H° and final 
states. 

Next we consider the process with one CP-even scalar and one CP-odd scalar ^4° in the final state, 
i- e - XX ~~ * h°A° and XX ~^ H°A°. The diagrams for these channels are shown in Fig. [8] and receive 
contributions from neutralino exchange in the t- and u-channel and from Z° and ^4° exchange in the 
s-channel. 

X? 71 A X? ^>^-:::-"' A x \.,l. <: -- 1 

X? ^ h,H xl^^ "--■h,H x\^^ '"•■-h,H 



Figure 8: Diagrams for the processes containing one pseudoscalar Higgs boson and one scalar Higgs 
boson in the final state. 

Finally, we consider the process XX ~* H + H~ . This proceeds by exchange of the two charginos in the 
t- and u-channel and by s-channel exchange of the Z°, hP and H° bosons, as shown in Fig. [9j 



Xl - 



Xi - 



h.H ..r-- H+ 



Figure 9: Diagrams for the process XX ~^ H + H . 
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4.2 Calculating cross sections at tree-level: MadGraph and CompHEP 

For all tree-level cross-section calculations we use both MadGraph/MadEvent [5] and CompHEP 

a- 

MadEvent is a multi-purpose, tree-level event generator which is powered by the matrix element 
creator MadGraph. It is a widely used tool in high energy physics to simulate events for collider ex- 
periments. Given a process and a specified model (SM, MSSM, 2HDM, etc), MadGraph automatically 
creates the amplitudes for all the relevant subprocesses and produces the mappings for the integration 
over the phase space. This information is passed to MadEvent, which calculates the cross section and 
may supply the user with unweighted events. These unweighted events are stored in the "Les Houches 
Event Files" (LHEF) format and may be passed directly to a Monte Carlo program (interface available 
for Pythia). However, for our purposes it is sufficient only to calculate the cross section. 
CompHEP is a package for automatic calculations of elementary particle decay and collision prop- 
erties at tree-level. The main idea prescribed in CompHEP is to make available passing on from 
the Lagrangian to the final distribution effectively with a high level of automation. CompHEP is a 
menu-driven system which makes it very user-friendly. It is divided into two parts: a symbolical and 
a numerical calculation. During the symbolical part, the user is able to specify and exclude diagrams 
in order to find out, how the different diagrams interfere. The numerical calculation may provide 
amongst others an energy dependent cross-section plot. 

4.3 Neutralino annihilation cross section at mSUGRA point SPSla 

In this section we study the cross section of the possible tree level neutralino annihilation processes at 
the mSUGRA benchmark point SPSla. This first study has mainly two purposes: We want to see (i) 
how the cross section behaves with the energy and (ii) which channels are dominant. Moreover, this 
study was used to cross check the results from both MadGraph and CompHEP in order to find out 
how they work and which assumptions and approximations are used for the calculation. 

4.3.1 mSUGRA benchmark point SPSla 

The mSUGRA parameters of the benchmark point SPSal, which is defined in |10| . are given in Table 
[TJ In Table [2] the mass values of some important particles are shown. Depending on the used SUSY 
spectrum calculator, the masses in MadGraph and CompHEP differ slightly but not significantly. The 
masses of the Higgs bosons A and H° are almost degenerate. 



m 


100 GeV 


m 12 


250 GeV 


tan (3 


10 


sign(^) 


+ 


^0 


-100 GeV 



Table 1: mSUGRA parameters of benchmark point SPSla. 
4.3.2 Energy dependence 

In the following, the energy dependence of the cross section of the most interesting tree level processes is 
presented. The <r(£ , )-plots are created by CompHEP in a zero-width-approximation of the exchanged 
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Particle 



Mass /GeV 




96.23 

176.97 

363.76 

381.73 

181.70 

379.94 

118.37 

394.23 

394.49 

402.53 



Table 2: Mass spectrum of mSUGRA benchmark point SPSla. 



Higgs bosons. Thus, it is possible to distinguish the resonance effects due to the s-channel exchange 
of ^4° and H , although they are almost mass-degenerate. Once we insert finite widths for these Higgs 
bosons, the distinguishablility of these resonances vanishes. 

The cross section has been calculated for random energy values with both CompHEP and MadGraph 
to do a cross check. The results are in reasonable agreement. 

Fermionic final states 

The energy dependence of the cross-section for the process XX ~* t + t~ in the range from 200 GeV 
(~ 2m x ) up to 1 TeV is shown in Fig. [TUl We see a clear resonance in the spectrum. A more detailed 
view on this resonance reveals two peaks at 394.2 GeV and 394.5 GeV. Comparing with the mass 
spectrum of SPSla in Table [2] we note that these peaks are due to the s-channel exchange of the Higgs 
bosons A and H°. In the region 2m x < < 2niA the cross section is of order (0.5 — 1.0) pb whereas, 
beyond the resonance, the cross section decreases as usual. 



Another important fermionic channel is the process XX ~> bb, whose cross section-energy dependence 
is shown in Fig. [TTJ At the resonance the cross section increases significantly up to order O(10 2 pb) 
and this process is therefore the dominant channel in this energy region. 




Figure 10: Cross section energy dependence 
of the process XX ~^ t + t~ . 



Figure 11: Cross section energy dependence 
of the process XX ~^ bb. 
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Final states containing gauge bosons 



Next we look at neutralino annihilation processes with gauge bosons in the final state, i.e. XX ~* 
W + W~ and XX ~^ Z°Z°. The cross sections of these processes are shown in Fig. [12] and [T3l 
respectively. Since there is no s-channel exchange of the pseudoscalar Higgs A there is only one 
resonance from the H° in the spectrum. In addition, there is a destructive interference between the 
s-channel diagrams of h° and H° which we study in detail in Appendix [5] for the process XX ~^ 

Z°Z°. 

The neutralino annihilation to gauge bosons is dominated by the process XX ~~ * W + W~. 




Figure 12: Energy dependence of the cross Figure 13: Energy dependence of the cross 

section of the process XX ~~ * W + W~. section of the process XX ~> 

Final states containing Higgs bosons 

Neutralino annihilation processes which contain Higgs bosons in the final state are only available 
at higher energies since the masses are much larger than m x . An exception may be the processes 
XX ~^ Z°h° and XX ~~ * hPhP, which contain the lightest Higgs in the final state. These channels 
are already open at lower energies and therefore contain resonances as well. The process XX ~> 

h°h° 

shows similar behavior as the process XX ~^ Z°Z° with a significant interference between the s-channel 
diagrams of h and H°. 

4.3.3 Comparing cross sections 

The cross sections for all tree-level neutralino annihilation processes have been calculated and are 
given in Table [3] for five different energies, i.e. 200 GeV, 360 GeV, 600 GeV, 1000 GeV and 10000 
GeV. The neutralino annihilation is dominated by processes with the leptons e, [i and r in the final 
state. The total cross section is about 1.7 pb at 200 GeV where the neutralinos are nearly at rest. At 
360 GeV where the energy is close to the mass of A and H° , the cross section increases up to around 

4.3 pb. In the resonance region, the neutralino annihilation is dominated by the process XX ~> bb. 

4.4 Parameter dependence 
4.4.1 Neutralino mass 

The dependence of the lightest neutralino mass m x on mo and mi/ 2 is sketched in Fig. [T4l and fT5| 
respectively, for tan/3 = 10, Aq = —100 GeV and fj, > 0. We see, that m x increases slightly with too 
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a[pb] 


Center-of-mass energy ^/s/GeV 


Final state 


200 


360 


600 


1000 


10000 


e + e~ 


4.89E-01 


8.79E-01 


4.72E-01 


2.04E-01 


2.37E-03 




4.89E-01 


8.79E-01 


4.72E-01 


2.04E-01 


2.37E-03 


T + T - 


5.44E-01 


9.57E-01 


4.82E-01 


2.06E-01 


2.38E-03 


v e v e 


3.14E-02 


6.69E-02 


4.03E-02 


1.86E-02 


2.30E-04 




3.14E-02 


6.69E-02 


4.03E-02 


1.86E-02 


2.30E-04 


V T V T 


3.18E-02 


6.73E-02 


4.05E-02 


1.84E-02 


2.30E-04 


uu 


8.13E-03 


4.32E-02 


6.10E-02 


5.27E-02 


1.35E-03 


dd 


6.41E-04 


2.08E-03 


3.46E-03 


3.28E-03 


9.10E-05 


cc 


8.13E-03 


4.32E-02 


6.10E-02 


5.27E-02 


1.35E-03 


ss 


6.41E-04 


2.08E-03 


3.46E-03 


3.28E-03 


9.10E-05 


tt 




5.56E-02 


9.93E-02 


7.49E-02 


1.45E-03 


bb 


6.76E-02 


1.13E+00 


3.11E-02 


5.64E-03 


1.05E-04 


w+w- 


9.19E-03 


2.54E-02 


1.71E-02 


1.15E-02 


5.46E-04 


zz 


3.12E-03 


1.01E-02 


5.02E-03 


3.36E-03 


2.44E-04 


Zh 




2.19E-02 


2.17E-02 


1.58E-02 


6.32E-04 


ZH 






2.66E-02 


1.40E-02 


1.94E-04 


ZA 






2.06E-02 


1.20E-02 


1.58E-04 


W + H~ 






2.37E-02 


1.28E-02 


2.25E-04 


W~H + 






2.37E-02 


1.28E-02 


2.25E-04 


hh 




2.88E-02 


4.91E-03 


2.52E-03 


2.11E-04 


HH 








5.69E-04 


2.50E-04 


hH 






3.06E-02 


1.37E-02 


1.05E-04 


H+H- 








4.45E-03 


6.83E-04 


Ah 






3.65E-02 


1.58E-02 


1.49E-04 


AH 








6.41E-03 


7.13E-04 


Otot 


1.71E+00 


4.28E+00 


2.02E+00 


0.99E+00 


1.66E-02 



Table 3: Cross section in pb of all tree-level neutralino annihilation channels at different energies 
for the mSUGRA-parameter point SPSal. 

whereas significantly with m^ 2 . The dependence on tan/3, Aq and sign(/x) is negligible. 




Figure 14: Neutralino mass dependence on Figure 15: Neutralino mass dependence on 

mo for mi/2 = 250 GeV, tan/3 = 10, Aq = "1-1/2 f° r m o = 100 GeV, tan/3 = 10, Aq = 

-100 GeV and n > 0. -100 GeV and fi > 0. 
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4.4.2 Total neutralino annihilation cross section 

In order to study the dependence of the neutralino annihilation cross section on the mSUGRA pa- 
rameter space scans have been performed in the mi^-mo-plane for fixed values of tan/3 = 20,50, 
Aq = -1000, 0, 1000 GeV and fi > for two different center-of-mass energies 500 GeV and 1000 GeV. 

For every parameter point, the particle spectrum was calculated by SoftSUSY [11] and decay widths 
were added by MSSMCalc [12J. Then, the cross-section calculation was performed by MadGraph. 
The m^-mo-plane was scanned in steps of 25 in each direction from 50 up to 1000. For the interpo- 
lation and plotting we used Root |13j . 

A parameter scan with incoming neutralinos nearly at rest would also have been interesting. However, 
MadGraph failed to calculate the cross section for a process with initial particles nearly at rest. But 
in fact, in the threshold regions, instead of looking at the cross sections one should consider the mean 
value of the annihilation cross section times the relative velocity, which is also the quantity to use in 
order to calculate the neutralino relic density. 

In Fig. [16] the total neutralino annihilation cross section at tree-level is shown for tan/3 = 20, Aq = 
and [i > and center-of-mass energies 500 GeV (left) and 1000 GeV (right). In addition, contour lines 
for neutralino masses of m x = 100 GeV, 200 GeV, 300 GeV and 400 GeV as well as for Higgs masses 
rriA ~ run = 500 GeV (left) and 1000 GeV (right) are included. The masses of the heavy neutral Higgs 
H° and of the pseudoscalar Higgs ^4° are almost degenerate nearly in the whole mSUGRA parameter 
space. Since the neutralino mass is not very sensitive to tan j3 and Aq, these contour lines are omitted 
in the following diagrams. 

For very small values of TOi/2> i.e. for light neutralino masses m x , the cross section is of order 0(1 pb) 
for y/s = 1000 GeV and of order O(10 pb) for y/s = 500 GeV. In this region the process XX ~* W + W~ 
is dominant, but there are also considerable contributions from XX ~^ Z°Z°. If too and TOi/ 2 are small 
annihilation processes with Higgs particles in the final state cannot be neglected, but with increase 
of too and m\i% these contributions decrease rapidly. If we approach a parameter region where the 
center-of-mass energy matches the mass of the Higgs particles ^4° and H°, the cross section increases 
significantly due to the resonance of these particles in the s-channel, which we discussed earlier in 
section [4.3.21 for the parameter point SPSla. Here, the process XX ~^ bb is dominant (rs 75 — 80%), 
but we get also a large contribution from XX ~> t + t~ (~ 16 — 18%) and considerable contributions 
from XX ->it (» 2%). 

Next we look at the parameter scan for tan (5 = 20, Aq = 1000 GeV and fi > for both energies 
500 GeV and 1000 GeV which are shown in Fig. [T71 Here, no specific increase of the cross section 
due to the process XX ~^ W + W~ appears for very low values of TOi/2- In fact, this process does not 
contribute significantly over all the scan region. In the resonance-region, the main contribution comes 
again from XX ~^ bb (« 80 — 82%), XX ~^ t + t~ (ss 15 — 17%) and XX ~^ & (~ 2%). For small values of 
too and TO1/2, SoftSUSY failed to generate a new parameter point due to theoretical inconsistencies. 
Thus, we denote this area as unphysical region. 

In Fig. [18] the total neutralino annihilation cross section at tree-level is shown for tan (3 = 20, 
Aq = —1000 GeV and \i > and center-of-mass energies 500 GeV (left) and 1000 GeV (right). 
Apart from the resonance, the main contributions come from processes with the leptonic final states 



11 



e + e~, /x + [i~ and r + r— and from the process XX ~> At the resonance we find the same behavior as 
before with a dominating bb production (~ 79 — 81%) and contributions from the channels XX ~^ t + t~ 
16 — 17%) and XX ~^ ti 2 — 2.5%). In the case y/s = 500 GeV the resonance region disappears 
into the unphysical region, i.e. the masses of ^4° and H° are larger than 500 GeV for all plotted points. 

Comparing the parameter scans in Fig. [TBI [T71 and [18] for the three values Aq = 0, 1000, -1000 
GeV, respectively, we may conclude that the neutralino annihilation cross section at tree-level is not 
very sensitive to the parameter Aq except for the slight shifting of the masses of A and H° and 
therewith of the resonance regions. 

Next, we look at the parameter plots for tan/3 = 50, \i > and A = 0, 1000, -1000 GeV in Fig. \W\ 
[20] and [21] respectively, for center-of-mass energies of 500 GeV (left) and 1 TeV (right). The masses 
of Aq and Hq decrease with tan /3, so that the resonance region shifts to higher values of tuq and m^/a- 
Therefore, the resonance region disappears from the scanned region for the plots with y/s = 1 TeV. 

For fixed parameters tan (3 = 50, Aq = and fx > (Fig. [T9|) the process XX ~* W + W~ con- 
tributes significantly for smaller than 100 GeV as seen before for tan/3 = 20, Aq = and fi > 
(Fig. [To]) . But in addition, there is also a large contribution from XX ~^ bb in this region which we did 
not see before in Fig. [T6j The resonance region for y/s = 500 GeV is again dominated by XX ~* bb 
(» 79 - 82%) and X X -> t + t- (» 18 - 20%). 

The total tree-level neutralino annihilation cross section for tan/3 = 50, Aq = 1000 GeV and fi > is 
shown in Fig. [20] At the resonance at center-of-mass energy y/s = 500 GeV the main contributions 
come from XX ~* bb (~ 80%) and XX ~^ t + t~ (~ 20%). Apart from the resonance, the most impor- 
tant channels are those with final states e + e~ , t + t~, bb and ti, the latter especially in the case 
of y/s = 1 TeV. 

In Fig. [2T]the parameter scans for tan/3 = 50, Aq = —1000 GeV and fi > for center-of-mass energies 
500 GeV and 1 TeV are shown. At the resonance, contributions come from XX ~* bb (~ 77 — 80%) 
and XX ~> t + t~ (~ 20 — 22%). In the remaining region, main contributions come from channels with 
final states e + e~, fi + fj>~, t + t~, bb and ti. 

Comparing the magnitute of the cross-section in Fig. [T9l [20l and [21] in the resonance region, we 
see that the cross section increases slightly with Aq. 

4.4.3 Conclusions of parameter dependence study 

We have seen that the resonances due to the s-channel exchange of ^4° and H° in the channels with 
massive fermionic final states lead to a significant increase of the total neutralino annihilation cross- 
section. In these resonance regions, the process XX ~* bb dominates (~ 75 — 82%). Other important 
contributions come from XX ~* t + t~ (~ 15 — 22%) and XX ~* (~ 2%). Since the masses of the 
Higgs bosons increase with itlq and mi2, channels with Higgs bosons in the final state only play a 
role for very low values of uiq and mm- Apart from the resonances, the most important channels are 
those with final states e + e~, r + r _ , bb and ti. 

The neutralino annhilation cross section seems not to be very sensitive to the parameter Aq. However, 
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Figure 16: Total tree-level neutralino annihilation cross section in my 2 -mo plane for tan/3 = 20, 
Aq = and \x > and center-of-mass energies 500 GeV (left) and 1000 GeV (right). Contour lines 
for neutralino masses of m x = 100 GeV, 200 GeV, 300 GeV and 400 GeV as well as for Higgs masses 
m A m H = 500 GeV (left) and 1000 GeV (right) were added. 



tree-level Vs=500 GeV mSUGRA -tanp=20 A o =1000 (i>0 log (u)/pb it^^ tree-level Vs=1000 GeV raSUGRA-tan|?=20 A o =1000 n>0 log(a)/pb 




100 200 300 400 500 600 100 200 300 400 500 600 700 800 900 1000 

m l/2 m l/2 

Figure 17: Total tree-level neutralino annihilation cross section in mift-Wo plane for tan/3 = 20, 
Aq = 1000 GeV and \i > and center-of-mass energies 500 GeV (left) and 1000 GeV (right). 
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Figure 18: Total tree-level neutralino annihilation cross section in m 1 / 2 - m -o plane for tan/3 = 20, 
Aq = —1000 GeV and ji > and center-of-mass energies 500 GeV (left) and 1000 GeV (right). 
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Figure 19: Total tree-level neutralino annihilation cross section in mirt-mo plane for tan/3 = 50, 
Aq = and \i > and center-of-mass energies 500 GeV (left) and 1000 GeV (right). 
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100 200 300 400 500 600 100 200 300 400 500 600 700 800 900 1000 

m l/2 m l/2 

Figure 20: Total tree-level neutralino annihilation cross section in mi^-mo plane for tan/3 = 50, 
Aq = 1000 GeV and fi > and center-of-mass energies 500 GeV (left) and 1000 GeV (right). 




Figure 21: Total tree-level neutralino annihilation cross section in my 2 -mo plane for tan/3 = 50, 
Aq = —1000 GeV and fi > and center-of-mass energies 500 GeV (left) and 1000 GeV (right). 
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for A = 0, there appears a region of increased cross section for very small m]/ 2 due to the process 
XX ~^ W + W~ . Since the masses of ^4° and H° depend on Aq there is a slight shifting of the resonance 
regions. 

The dependence on tan (5 seems to be more delicate, since the masses of the Higgs bosons A and 
H° are more sensitive to this parameter. Thus, for a higher tan/3, the resonance region appears for a 
given center-of-mass energy at higher values of mo and m-i/2- 

We want to note that, for cosmology, those parameter regions where 2m x ~ uia are of greater impor- 
tance. Large regions of the mSUGRA parameter space predict much larger values for the neutralino 
relic density than the values given by WMAP. To produce the observed dark matter relic density 
an enhanced neutralino annihilation is needed. One possibility would be an accidental degeneracy 
2m x ~ niA- Assuming CP-invariance, the s-channel exchange of the A boson would still be open at 
threshold, while the exchange of the CP-even H° would be suppressed. 

5 Summary and outlook 

We studied the cross section of all tree-level neutralino annihilation processes concerning energy de- 
pendence and channel dominance for the mSUGRA benchmark point SPSla. In processes which 
contain contributions from s-channel exchange of a neutral heavy Higgs H° and/or a pseudoscalar 
Higgs boson A , these diagrams lead to resonances when the exchanged Higgs bosons are produced 
on-shell. Since the Higgs bosons couple to the mass, those channels with third generation fermions in 
the final state, particularly the process XX ~> bb, dominate the total cross section in this resonance 
region. Apart from the resonance region, the neutralino annihilation cross section receives the largest 
contributions from processes with e + e~ , fjr and r + r _ final states. 

We studied the mSUGRA parameter dependence of the neutralino annihilation cross section by scan- 
ning the m 1 /2- m o plane for fixed values tan/3 = 20, 50, Aq = —1000, 0, 1000 GeV and fi > for 
two different center-of-mass energies -y/s = 500 GeV, 1000 GeV. The parameter scans have shown 
that the dominance of the process \X ~^ bb in the resonance region is not sensitive to the mSUGRA 
parameters. The behavior apart from the resonance, i.e. that the main contributions come from 
channels XX ~^ e + e~, fi + fi~, t + t~, does not vary significantly with the mSUGRA parameters either. 
However, there appear some regions in the parameter space where the processes XX ~> W + W~ and 
XX ~* ti provide considerable contributions. For Aq = 0, there is even a small region at low values of 
rrii/2, where the process XX ~^ W + W~ is dominating. Since the masses of the Higgs bosons ^4° and 
H° depend on the mSUGRA parameter, the resonance region shifts in the mo-m^-phine for a given 
center-of-mass energy. 

After this study on neutralino annihilation processes at tree-level we started to investigate higher- 
order corrections by calculating the cross section up to one-loop-level with FormCalc |9j. We also 
studied neutralino annihilation to the two photon final state where the leading order happens to be at 
loop-level. The tree-level results, obtained by FormCalc, are in good agreement with those derived 
with CompHEP. At loop-level, the overall form of the cross section plots seems reasonable. How- 
ever, since counter-terms have not yet been included for the MSSM in FormCalc, we are not totally 
certain about our results and thus decided not to present them here. 



16 



Appendix A 

Detailed view on the Feynman diagrams contributing to 

In this appendix, we study the contributions and interferences of the six different Feynman diagrams 
for the process XX ~* 

Z°Z° given in Fig. [22] at the mSUGRA point SPSla, as an example. For this, we 




Figure 22: Six Feynman diagrams are contributing to the process XX ~^ Z°Z°: exchange of all 4 
neutralinos in the t-channel (n= 1,2, 3,4) and s-channel exchange of the neutral Higgs bosons h and 
H°. 

use the possibility to exclude specific diagrams in CompHEP in order to extract the contributions of 
those diagrams we are interested in. Then, we plot the cross section energy dependence in the region 
from 200 GeV up to 1 TeV. In general, the cross section is proportional to the scattering amplitude 
squared, i.e. the sum of all diagrams squared: 

a oc \M\ = Vy. diagrams (4) 

Hence, we have all in all 21 terms contributing, where 6 are squared diagrams and 15 are interference 
terms. In the following, we show the most important contributions. 



In Fig. [23 we show again the whole cross-section of the process XX ~^ Z°Z° with all diagrams 
contributing. As mentioned before, there is a resonance at 394.5 GeV due to the s-channel exchange 
of H°. In addition, there happens to be a destructive interference, we want to find out about. 




Figure 23: Cross section energy dependence of the process XX ~* Z Z with all diagrams contributing. 

We start by looking at the six squared single diagram contributions. In Fig. Othe cross section energy 
dependence of one single t-channel diagram is shown. Since the slope is for all four t-channel diagrams 
similar, we omit here to show all four plots. But in contrary, the magnitude of these diagrams differs 
significantly. An overview of the magnitudes is given in Table HI Note, that these values are only 
given in order to get a feeling for the order of magnitudes. However, we see that the main t-channel 
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200 



Figure 24: Contribution of single t-channel diagram \X ~^ (Xn) ~^ Z°Z°. The slope is for all four 
neutralinos similar while the magnitude differs significantly (see text). 





exchanged neutralino 


Energy 


x [ l 


xl 


xi 


xi 


200 


3E-06 


3E-06 


2E-04 


1E-05 


300 


7E-06 


4E-05 


2E-03 


1E-04 


1000 


1E-05 


2E-04 


2E-02 


3E-03 



Table 4: Cross section in pb for the four single t-channel diagram contributions with neutralino 
exchange for three different energy values yfs = 200 GeV, 300 GeV, 1 TeV. 

diagram contributing is via the exchange of x-j- 

The cross section contribution of the single s-channel diagrams with exchanged h° and H° are shown 
in Fig. [25] and [26] respectively. While the cross section contribution of the diagram with /i°-exchange 
increases monotonically, the ff°-exchange diagram contribution has a broad resonance and decreases 
for higher energies. 




Figure 25: Contribution of the s-channel di- Figure 26: Contribution of the s-channel di- 

agram XX -> (h°) -> Z°Z°. agram X X -> (H°) -» Z°Z°. 

Now, we look at the 15 interference terms, that we separate into t-channel interferences, s-channel 
interference and s-t-channel interferences. Note, that although we discuss in the following the inter- 
ference terms, the squared diagrams are also included in the plots. We select (two) specific diagrams 
and calculate the cross-section from these diagrams, which of course also includes the squared single 
diagram contributions. 

The t-channel interference, i.e. the six interference terms between the t-channel diagrams exchanging 
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neutralinos, is reasonably simple. All interferences are dominated by the diagram exchanging x% if 
it is involved. The slopes of the t-channel interferences are similar to those of the single t-channel 
diagram contributions in Fig. 1241 The t-channel interference terms are almost negligible compared to 
the squared term of the diagram exchanging x®- This can be seen in Fig. [271 where the total t-channel 
contribution, i.e. all four diagrams exchanging neutralinos, is shown. 




Figure 27: Contribution from all four t- 
channel diagrams XX ~^ (Xn) ~^ Z°Z°. 



Figure 28: Interference of the s-channel di- 
agrams XX ~ > (h ) — * Z°Z° and XX ~~ > 
(H°) Z°Z°. 



The interference between the s-channel diagrams with exchange of h° and H° is shown in Fig. [28l 
Here, we see a significant interference at about 400 GeV, i.e. near the resonance of the heavy neutral 
Higgs boson H°. This gives a first hint, where the large interference in the cross section of the process 
xx -» z°z° comes from. At higher energies, the s-channel interference contribution seems to be 
constant or increasing slowly. Thus, the decrease of the cross section at higher energies results from 
interferences between t- and s-channel. 



In Fig. [29] - [32] we study the interference between the s-channel diagram exchanging a light Higgs 
boson h° and one of the four neutralinos exchanging t-channel diagrams. The cross section for those 
interference terms involving t-channel diagrams with x1, X2 an d X% are quite similar in slope and mag- 
nitude, while the cross section for the interference term involving X3 is about one order of magnitude 
lower and features a clear maximum at about 270 GeV and a minimum at about 650 GeV. 




Figure 29: Interference of the s-channel dia- 
gram XX ~^ (h°) — ► Z°Z° and the t-channel 
diagram XX -> (x?) -> Z°Z°. 



Figure 30: Interference of the s-channel dia- 
gram XX ~^ {h°) Z°Z° and the t-channel 
diagram XX -> (x°) Z°Z°. 
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Figure 31: Interference of the s-channel dia- 
gram XX ~^ (h°) — ► Z°Z° and the t-channel 
diagram XX -> 



Figure 32: Interference of the s-channel dia- 
gram XX ~^ (" ) — ► Z°Z° and the t-channel 
diagram XX ~^ ~> -Z^ - 



Finally, we look at the interference term contributions between the s-channel diagram exchanging a 
heavy neutral Higgs boson H° and one of the four neutralinos exchanging t-channel diagrams, shown 
in Fig. [33] - [361 In all plots we see a resonance at around 400 GeV, i.e. at the mass of the Higgs 
boson H°, In the interferences containing x?, X2 anc ^ X4 i n Fig- l33 | 1341 and I36 [ respectively, we have 
first the maximum and then a minimum (or a monotonically decreasing slope in the case of Fig. I331) . 
However, for the interference between the diagrams exchanging H° and X3 we have the opposite, first 
a minimum and then a maximum. 




Figure 33: Interference of the s-channel dia- 
gram XX ~^ — ► Z°Z° and the t-channel 
diagram XX -> (x?) ^ Z°Z°. 



Figure 34: Interference of the s-channel dia- 
gram XX ~^ — ► Z°Z° and the t-channel 
diagram XX (x°) ^ 



Concluding this section, we may say that the main contribution to the large interference in the 
cross section of the process XX ~^ Z°Z° is due to the interference between the s-channel diagrams 
exchanging h° and H°. The interference between s-channel diagram containing H° and the t-channel 
diagrams with neutralino exchange may also have some impact on this behavior in the resonance 
region. Moreover, the t-channel diagram containing ^3 differs in its behavior from the other t-channel 
diagrams. However, it is difficult to understand all these features only by studying the cross sections of 
a selected subset of diagrams. Thus, finding out the main contributing diagrams in a process contain- 
ing overall six diagrams is not at all a trivial task. In principle, in quantum mechanics, all diagrams 
interfere. Any subdivision (as here) is artificial. Still it gives good insights to understand this complex 
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Figure 35: Interference of the s-channel dia- 
gram XX ~^ — ► Z°Z° and the t-channel 
diagram XX (x°) -> 



Figure 36: Interference of the s-channel dia- 
gram XX ~^ — ► and the t-channel 
diagram XX (x°) -> 



phenomenon. 
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